Using a qubit to measure photon number statistics of a driven, thermal oscillator 
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We demonstrate theoretically how photon number statistics of a driven, damped oscillator at 
finite temperature can be extracted by measuring the dephasing spectrum of a two-level system 
dispersively coupled to the oscillator; we thus extend the work of Dykman [lj and Gambetta et al. 
0| • We carefully consider the fidelity of this scheme- to what extent does the measurement reflect the 
initial number statistics of the mode? We also derive analytic results for the dephasing of a qubit by a 
driven, thermal mode, and compare results obtained at different levels of approximation. Our results 
have relevance both to experiments in circuit cavity QED, as well as to nano-electromechanical 
systems. 
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I. INTRODUCTION 

There has been a recent explosion of interest in the 
properties of quantum oscillators. On one hand, nano- 
electromechanical systems (NEMS) (systems comprising 
of a micron scale oscillator coupled to an electronic con- 
ductor) have been used for near quantum-limited position 
detection, and are thought to be good candidates for the 
preparation of exotic quantum states and for use in quan- 
tum control experiments On the other hand, there 
has been considerable success in studying circuit QED 
systems, in which a superconducting qubit is strongly 
coupled to a single mode in an electromagnetic resonator 
[1, IH @ • Experiments in these systems have been able to 
enter the strong coupling regime of cavity QED, achiev- 
ing such milestones as the resolution of the vacuum Rabi 
splitting Q. 

A clear goal in both these classes of systems is to reach 
a regime where the quantum nature of the resonator is 
unambiguously manifest in measured experimental quan- 
tities. An example of this would be the observation of 
features associated with the discrete energy states of the 
oscillator. This was recently achieved in spectacular fash- 
ion in the circuit QED system with a so-called "number 
splitting" experiment 0]. It was experimentally shown 
that for a strong dispersive coupling between the qubit 
and the mode, the qubit dephasing spectrum (i.e. the 
time-dependence of the qubit 's off-diagonal density ma- 
trix element) could be used to probe the photon number 
distribution in the cavity mode. Theoretically, this effect 
was first described by Dykman for the case of a undriven, 
thermal mode More recently, Gambetta et al. 
described the effect for the case of a zero-temperature 
driven mode, a case more relevant to the experiment of 
Ref. 0. 

In the case of NEMS systems, there have been sev- 
eral theoretical proposals for measuring discrete number 
states, but to date these have been difficult to implement 
experimentally [1, [lj| • Given the success of the circuit 
QED number splitting experiment Q , it is natural to ask 
whether a similar approach could be taken with NEMS 
systems. In principle, the answer is yes. The coupling 



of nanomechanical oscillators to superconducting qubits 
has been discussed extensively [ill EH EH an d is being 
actively pursued experimentally. One could thus easily 
envisage a setup where a superconducting qubit is cou- 
pled to both a electromagnetic cavity and a nanomechan- 
ical oscillator 0, [T3 ]. The mechanical oscillator would 
be strongly coupled to the qubit and would generate the 
number splitting effect, while the electromagnetic cav- 
ity would be weakly coupled and would be used simply 
to read out the qubit dephasing spectrum. The idea of 
using the approach of Ref. [7| to detect number states in 
a nano-mechanical oscillator is highly attractive. From 
an experimental perspective, it is not as demanding as 
previous proposals (e.g. strong non-linearities in the os- 
cillator are not required). In addition, we will see that 
such an experiment offers the possibility of providing an 
unambiguous signature of the discrete number states of 
the oscillator; no extensive post-experiment interpreta- 
tion is required. 

The above discussion notwithstanding, there remain 
unresolved theoretical issues with the number splitting 
effect that should be addressed before contemplating a 
NEMS version of the experiment. These issues are also 
relevant to the circuit-QED experiment of Ref. 0. First, 
as NEMS oscillators do not typically have frequencies 
much larger than temperature, one should now consider 
both the effect of finite temperature and a driving force; 
this regime is not addressed by existing theory [l|, Q . In 
this paper, we calculate analytically the qubit dephasing 
spectrum for a qubit dispersively coupled to a driven, fi- 
nite temperature oscillator, using a conceptually simple 
Wigner function approach. We recover the previous re- 
sults in the appropriate limits. Our results also provide 
analytic expression for the long-time oscillator-induced 
dephasing rate of the qubit in the driven, thermal case 
(c.f Eqs. dU]) and (|tj5|) ). Note that the dephasing rate in 
the absence of a drive was considered recently by Serban 
et al. in Ref. GJ. 

A second important issue we address here is the ques- 
tion of measurement fidelity: to what extent is the qubit 
dephasing spectrum a faithful measurement of the mode's 
number statistics? It was already pointed out in Refs.0,0 
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that for a generic driving force, there was not a simple 
connection between the observed peaks in the qubit de- 
phasing spectrum and the photon number distribution in 
the mode. We demonstrate that in general, back-action 
effects arise from the finite damping of the oscillator; 
these imply that the dephasing spectrum of the qubit 
will fail to perfectly reflect the initial photon distribu- 
tion of the oscillator. Nonetheless, we demonstrate that 
for a sufficiently strong oscillator-qubit coupling, one can 
reliably extract the initial photon distribution of the os- 
cillator. 



a scheme, one attempts to excite the qubit from its 
ground state via a time-dependent field at frequency u> r f 
which couples to a x . Within Fermi's golden rule and 
the quantum regression theorem [l7| , the absorption rate 
T abs (aj rf ) is given by: 



T abs (cj rf ) = \A\ 2 dte+to't* (a-(t)tr+(0)) 



2|A| 2 Re 



dte +iu) rft 



uo 



(6) 



II. HEURISTIC DISCUSSION OF THE 
NUMBER SPLITTING EFFECT 

We consider a qubit which is dispersively coupled to a 
harmonic oscillator. Ignoring for the moment the damp- 
ing and driving of the oscillator, the Hamiltonian of the 
system is given by: 

H = H + H in t (1) 
H = ^L az +hn^a+^j (2) 

H int = XHil (a^a + a z (3) 

Here, Ml q b is the energy splitting of the qubit, f2 is the 
oscillator frequency, and A is a dimensionless measure of 
the interaction strength. As discussed, e.g. in Ref. 0, an 
interaction of this type can arise from a Jaynes-Cumming 
type coupling in the limit where the frequency differ- 
ence \fl q b — f2| is large (i.e. one transforms the starting 
Hamiltonian to eliminate the first-order qubit-mode in- 
teraction; the dispersive interaction above then results 
at second order in the original coupling strength). It has 
recently been suggested [15| that it is important to keep 
non-rotating terms in the interaction Hamiltonian Hi nt 
(i.e. terms proportional to a 2 and (a ! ) 2 ). We find that 
such terms play no role in the weak damping, weak cou- 
pling limit of interest (i.e. A <C 1,7/T2 <C 1); we discuss 
this further in what follows. Note that in this limit, the 
effects of the qubit-oscillator interaction may still be sig- 
nificant, as Afi/7 can be large. 

The quantity we wish to study is the qubit 's off- 
diagonal density matrix element. Letting | T)> I I) denote 
o z eigenstates and p the density matrix of the oscillator- 
qubit system, we define: 

p n (t) = p n (t) = (t \p(t)\ I) (4) 
Pn(t) = Tr osc p n (t) (5) 

Note that Pn(t) is indeed a measurable quantity; one 
could imagine extracting it using state tomography, 
as has recently been demonstrated for superconduct- 
ing qubits [l6| . An alternate method, as discussed in 
Refs. 0, Q and used in the experiment of Ref. 0, is to 
measure pn(t) via an absorption experiment. In such 



Here, the second line assumes the quantum regression 
theorem, and A is a matrix element proportional to the 
strength of the driving field. 



A. Number splitting interpretation of ptx(t) 

The simplest expectation for the behaviour of p-\i(t) 
follows from the observation that a) this quantity should 
oscillate at the qubit splitting frequency, and b) the dis- 
persive coupling in Eq. ([3]) implies qubit splitting fre- 
quency depends on the number of quanta in the mode. 
This leads to the expected form: 

oo 

PU (*) = e~ mqbt E P{n)e- 2lXn ^y Pn (0) (7) 

71=0 

where P(n) is the number distribution of the harmonic 
mode. Eq. ([Jj would of course be exact in the absence of 
damping or forcing of the oscillator. It tells us that the 
time- dependence of p-\i at a fixed coupling strength con- 
tains information about the oscillator photon number dis- 
tribution at t = 0: if one Fourier-transforms p-fj.(i), one 
expects a set of discrete peaks with a spacing Auj = 2AO, 
with the weight of each peak yielding the probability 
P(n) of the corresponding number state. 

In the more general case where we have to contend 
with damping and forcing of the oscillator, Eq. ([7]) is 
no longer a priori valid. The interaction Hi nt now fails 
to commute with the total oscillator Hamiltonian, and 
thus there will be a back-action associated with the mea- 
surement. Heuristically, turning on the coupling to the 
qubit can now itself cause a change in the number statis- 
tics of the oscillator. For example, the coupling to the 
qubit causes the oscillator's frequency to change; at non- 
zero damping, this will cause it to re-equilibrate with the 
bath and thus change its number distribution. We can 
nonetheless attempt to use Eq. ([7]) to extract the number 
statistics of the oscillator (i.e. use the weights of peaks in 
the spectrum of p^i to extract the distribution of number 
states in the oscillator). We will thus refer to Eq. ([7]) as 
the "number splitting" interpretation of pn; in what fol- 
lows, we will assess the validity of such an interpretation. 
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Full counting statistics interpretation of pn(X) trinsic dissipation of the oscillator; we choose the den- 
sity of states of these modes and the couplings Vj to 



A second, somewhat different way to connect pn(t) to 
number statistics is to use the coupling dependence of this 
quantity (as opposed to its time dependence) to extract 
the oscillator number statistics. This is the procedure 
used in defining the full counting statistics of transmit- 
ted charge in mesoscopic electron conductors 18]. The 
application of this approach to the present system will 
be discussed in a forthcoming work [2(| ■ 



yield a standard Ohmic dissipation. Note we are con- 
sidering an ideal case where the qubit has no intrin- 
sic (i.e. oscillator-independent) relaxation or dephasing. 
Such intrinsic processes can be accounted for by multiply- 
ing p-\i(t) for the ideal case by an extra decaying factor 



exp 



i 

2T, 
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III. EQUATION OF MOTION FOR p n 



We now wish to calculate p^i (t) in the presence of both 
oscillator damping and driving. The Hamiltonian of our 
system H thus takes the form: 

H = H + H int -mF(t)(a + a^) (8) 

+ (a + at) V J ( b i + b l) + E ^ ( b h + \ 

j j 

F(t) is the external drive on the oscillator. The bj modes 
correspond to a thermal bath responsible for the in- 



We are interested in the evolution of Pn(t), which is 
an operator in the oscillator's Hilbert space. We are also 
interested in the weak damping limit 7 <C f2 which is rel- 
evant to experiment. In this limit, one can start with an 
exact path-integral description and then rigorously de- 
rive a local-in-time equation of motion for p^; this is 
in complete analogy to the derivation of the usual weak- 
coupling Caledeira-Leggett master equation pjj . This 
equation can also be rigorously justified in the high tem- 
perature ksT 3> Ml limit. Setting Sl q b = (as it plays 
no role in what follows), one has: 



at 



PU 



:[H sc,Pn] 



F(t) 



2h 



%,{Pi Pn}} 



D 



[x, [x,pu]\ - tH h osc, pu} 



(9) 



r 



Here, H osc — p 2 /(2m) + m£l 2 x 2 /2, 7 is the damping rate 
of the oscillator due to the bath, and (/3 = l/(ksT)) 



D = mj^- coth(/3/ifi/2) 



(10) 



is the corresponding momentum diffusion constant. We 
also use Axq and Apo to denote the zero-point variances 
of the oscillator. Note that unlike the standard quantum 
optics master equation (used in Refs. [l], 0), the above 



equation has counter-rotating terms (i.e. proportional to 
a 2 and (a^) 2 ) and is hence not in Linblad form. We will 
comment more on this difference in what follows. 

To proceed, we will represent p-^ by its Wigner func- 
tion: 



W(x,p; t) 



1 

2vrft 



dye-^/Kpuix + ^x 



y 



t) (11) 



Writing Eq. ([9]) in terms of the Wigner function, one finds 



d t W(x,p;t) = 



P a 
— o x 
m 



mQ 2 x 



Ml 

Ax ' 



F(t) \d p + 1 d p -p + Ddp 



?L+mn 2 x 2 -%( -81 
m 4 



W{x,p- 1) 
(12) 



r 



The A-independent terms above are the usual terms in a 
classical Fokkcr-Planck equation for a damped, driven os- 
cillator, while the terms proportional to A yield distinctly 
non-classical contributions. As the additional quantum 
terms are at most quadratic in x,p and their derivatives, 



we may solve this equation in the same way one would 
solve the classical Fokkcr-Planck equation: one makes a 
Gaussian ansatz. Note that a similar approach was used 
recently by Serban et al. [l5[ to study the dephasing of 
a qubit by a thermal, damped mode. 
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Thus, we write: 

W(x,p) = ^-e ik *e^W(k,q) 
and make the Gaussian ansatz: 



W(k,q;t) 
W(0,0;0) 



gii/ gikAxo-x+iqApo-p x 



(13) 



(14) 



exp 



-{k 2 {Ax f<x x + q 2 (Ap ) 2 a p 



-2kq(Ax )(Ap Q )a : 



xp) 



The Gaussian is parameterized by six time-dependent 
parameters: the means x(t), p(t), the variances a x (t) 
and cr p (t), the cross-correlation a xp (t) and the overall 
phase v(t); note that we have scaled all quantities by 
the zero-point variances of the oscillator. Substituting 
the above ansatz into the appropriately transformed ver- 
sion of Eq. ([9]), we find that it is satisfies the equation 
as long as the Gaussian parameters satisfy the following 
equations (we set 17 = 1 from this point onwards): 



dx 

~dt 
dp 

~dl 

da x 

~dt 

d<j p 

~dT 

d(j X p 

dt 
dv 

~dt 



p~i\ (a x x + a xp p) 

—x + 2F(t) — jp — iX (<jpp + a xp x) 

-2a xp - 2^(a p 



. , p -i\ {a 2 x + alp - l) 



f) - i\ (a 2 



xp 



i\a xp (a x + a p ) 



A 



P 



Here, we have defined: 

f = coth/3/ift/2 



2n e 



1 



(15) 
(16) 
(17) 
1X18) 
(19) 
(20) 

(21) 



The above equations have a simple form. The A- 
independent terms describe the usual dynamics of a 
driven dissipative oscillator, and have a completely classi- 
cal form. The A-dependent terms act as quantum source 
terms, and cause the means and variances of the Gaus- 
sian to acquire imaginary parts. From Eq. (|20|) . we find 
that the decay of p^ is determined from these quantities: 



PTl(O) 



with 



Ai(i) = exp 
A 2 {t) = exp 



W(k =0,q = 0-t) 
W(k = 0,q= 0;0) 
Ai(i)-A 2 (t) (22) 



-4 / dt'(a x (t')+a p (t'))} (23) 
z Jo 

-i±J*dti(x 2 (t>)+p 2 (t>))\ (24) 



The two factors above correspond to two distinct con- 
tributions to the decay of p^, and hence to the num- 
ber splitting effect. The factor Ai(t) describes a purely 
thermal effect: from Eqs. p7|) and (fT5)) . we see that 
a x and a p acquire negative imaginary parts, causing 
Ai(i) to decay from its initial value of 1. This effect 
only occurs at finite temperature, as at zero temperature 
"x(0) = Cy(0) = T = 1 and are hence time independent. 
Note that Ai(t) is completely independent of the means 
of the Gaussian (i.e. x and p) and hence of the driving 
force F{t). 

In contrast, the factor A 2 (i) in Eq. (|2"2"|) describes the 
effect of a finite drive on p^: if the means x and p are 
non-zero, Eq. (|15|l - (fTB|) indicate that they will acquire an 
imaginary part and cause this factor to decay. 



IV. DEPHASING SPECTRUM IN THE PURELY 
THERMAL CASE 

We first consider the case of zero drive (i.e. F(t) = 0), 
and take the oscillator to initially be in a thermal state: 
x, — p — 0, a x = Up — T, and a xp = 0. Note that this 
corresponds to an oscillator in equilibrium with frequency 
O; one might argue that for the absorption experiment, 
the initial condition should correspond to an equilibrium 
oscillator with frequency 17 — A. This difference is easy 
to include; in the weak coupling limit of interest (A <C 
1), the corrections are lower-order in A and hence not- 
significant. 

As there is no drive, Eqs. (|15[) and p^]) imply x = 0, 
p = for all time, and thus the factor A 2 (t) = 1. The 
evolution of pn is the completely determined by the fac- 
tor Ai(t). To proceed with the analysis, it is convenient 
to write: 



<j x = 1 + s + As 
a v = 1 + s — As 



(25) 
(26) 



The equations for the Gaussian parameters thus become: 

(27) 



J = -zA[2s + s 2 + (As; 



2 + a lp] 



(28) 



dAs 

~~dT 

do X p 
dt 



—7(5 - As) + 2"fn eq 
-2iA(l + s)(As) 

+ 7 (s - As) - 2jn eq + 2a xv (29) 
— 7<7 xp — 2As — 2iXa xp [1 + s] (30) 



As we start in a thermal state, both As(0) = and 
&x P {0) — 0; further, if we took the limit 7^0, these 
quantities would remain zero for all time. Thus, in the 
weak damping (7 <C 1) limit of interest, it will be an 
excellent approximation to set As = a xp = in Eq. 
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this results in the approximate equation: 

-iX [2s + s 2 ] -7(s- 2n eq ) (31) 



ds 
dt 



Within this approximation, A\(t) is completely deter- 
mined by Eq. (f3Tj) and Eq. (|27|) . Note crucially that our 
approximation does not amount to doing perturbation 
theory in A or 7, as terms to all orders in the ratio A/7 
are retained. Also note that if one had started by using 
the quantum optics version of Eq. ([9]) (i.e. drop all a 2 and 
(a^) 2 terms), one would immediately find Eq. pip with- 
out any additional approximation. We thus see that the 
differences between the Linblad and non-Linblad forms 
of the master equation play no role in the weak coupling, 
small damping limit of interest. In what follows, we will 
explicitly discuss the correction terms that arise when 
one uses the full equation Eq. (|28j) versus the approxi- 
mate version Eq. (|31[) . 



A. Results for dephasing spectrum 



We can now solve directly Eq. (f3"Tj) for s. We find: 



, . H 1 + Me 



-2i\at 



2A 



I -Me 



-1i\at 



where 



M 



2e-*> + (e-/*> _ 1) (a-(l-&)) 



2 + (e-^ - 1) (1 - 



(32) 



(33) 



(34) 



2A/ 



Note that at T = 0, we have simply a = 1 — 27/ (2A) and 
M = 0. 

Combining Eqs. and ||32J) we find: 



m(o) 



= Ai(t) 



1 - M 



(35) 



1 - Me- 2iAat 

OC 

= e~< t/2 (l- M)J2 Mne ~ i(2n+1)Xat ( 36 ) 

Note that the result above coincides exactly with Eq. (11) 
in Ref. U. Unlike that work, we did not start by making 
a rotating-wave approximation in the bath-oscillator in- 
teraction, but instead took a weak-damping limit. 

Consider first the 7 — > limit. In this limit, one finds 
M = e _/3n and a = 1. Eq. J36]) thus yields p n (t) in 
the form of Eq. {3 with a thermal number distribution. 
As expected, is a completely faithful measure of the 
oscillator number statistics in this limit, as there is no 
back-action: at zero damping, the mode has no way to 
re-equilibrate its number distribution in response to a 
change in its frequency. 



In the case of small but finite 7, we see that \M\ will 
play the role of the Boltzman weight e - ' 3 ", and that a 
acts to renormalize the coupling A. For small damping, 
the Fourier spectrum of pn(t) will again consist of a set 
of equally-spaced discrete peaks occurring at frequencies 
uj n . The peaks will now have a Lorentzian form, each 
having a non-zero width j n . To leading order in 7 we 
have: 



Lu n = (2n + 1)A 



1+ \ 









(37) 
(38) 



7n = 7 i n eq + n(l + 2n eq )} 

Further, to order j 2 , the weighting factor \M\ is given 
by: 



\M\ ~ e 



-pn 



(39) 



We can use \M\ to define an effective temperature T e ff 
via \M\ = e ^ fi °/( fc '> T =//). This is the temperature one 
would associate with the oscillator based on using Eq. Q 
to interpret the peaks in the spectrum of pi±{t). To lead- 
ing order in 7/A, we have: 



T. 



eff 



= T x 1 



7 



4A 2 hVL 



(1 + 2n eq ) 



(40) 



Thus, in the weak-but finite damping limit, the proba- 
bility distribution extracted from the time dependence of 
Pll at fixed coupling strength A ^> 7 (i.e. using Eq. (JT])) 
is very close to the initial, equilibrium photon number 
distribution of the oscillator. Note that there is a one-to- 
one correspondence between the peaks in the spectrum 
and the peaks expected from Eq. {?]). The main effect 
of the damping is to give the peaks associated with each 
number state a finite width in frequency. There is also 
deviation between the weight of each peak and a pure 
thermal weight; the finite damping makes the oscillator 
appear somewhat colder than it should. On a heuristic 
level, it is not difficult to understand why we have a cor- 
respondence. The time-scale for back-action effects to 
set in is ~ I/7- this is how long the oscillator takes to 
re-equilibrate its number distribution given the change 
in frequency brought on by the coupling to the qubit. 
For A > 7, the qubit effectively "measures" the number 
statistics well before they have a chance to be changed 
by the back-action. 



B. Results for oscillator-induced dephasing rate 

It is also interesting to consider the limiting behaviour 
of p-\i(t) in the long-time limit, t 3> I/7. We can define 
the long-time dephasing rate T Vtt h and frequency shift 
A t h of the qubit via: 



Hh 



(p,th 



= — Im lim 

t — >OG 

= —Re lim 

t— >oo 



lo gPTiO) 
t 

logp u (i) 



(41) 
(42) 
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This yields: 



A 



Hi 



■ ip,th 




(43) 



(44) 

We find that the thermal oscillator-induced dephasing 
rate may be written as 7 times a function which depends 
only on the ratio A/7 and the temperature; in typical 
experimental situations, one has both 7^1 and A <C 
1, but the ratio A/7 can be arbitrary. Note also that 
the thermal dephasing rate of the qubit vanishes at zero 
temperature. Expanding in 7/ A we have: 



This confirms our prior claim that the error introduced 
by the approximation of Eq. (|31[) is small in the weak 
damping 7 <§; 1 limit. 

Finally, we have also calculated the thermal dephasing 
rate for a qubit-oscillator interaction Hamiltonian which 
contains counter-rotating wave terms, i.e. : 



Hi, 



(50) 



The thermal dephasing rate of the qubit for the inter- 
action Hamiltonian Hi n t, 2 was previously considered by 
Serban et al. [la |. Using a completely analogous pro- 
cedure as the above, we can find the oscillator induced 
dephasing rate of the qubit with this modified interac- 
tion; results are provided in the appendix. As expected, 
we find that corrections to the result of Eq. (|44|) are neg- 
ligible in the 7, A <C 1 limit of interest. 



V. DEPHASING SPECTRUM IN THE 
PRESENCE OF AN OSCILLATOR DRIVE 



At* 



A+^n e(J (l + n e9 ) + 0(7 4 ) 

7 3 - - 11 

in eq ~ -^2 n eq (1 + n eq ) (1 + 2n eq ) + O [ — 



(45) 

,5- 

(46) 



In contrast, in the high-temperature limit n eq 3> 1, 
find for the dephasing rate: 



r 



ip,th 



V7^ 



(47) 



Finally, expanding T v-t h as a function of the coupling 
yields the usual Fermi's golden rule result: 



F 



ip,th 



A 



2 4n eg (l 



7 



X 2 S n (u; = 0) 



(48) 



where S n {w — 0) is the zero-frequency charge noise of 
the oscillator at zero coupling (i.e. A = 0). Note that the 
high-tempcraturc limit of the Golden rule result and the 
full answer are not the same! This can be understood 
from the fact that higher order terms in A become in- 
creasingly important at higher temperatures; the Golden 
rule result is only valid if Xn eq /j <C 1 is satisfied. 

Eqs. (|44| - (l47p for the thermal dephasing rate are 
valid within the small-damping approximation made in 
Eq. pip , or equivalently, within a rotating- wave descrip- 
tion of the bath-oscillator coupling. One can also find 
^ip.th without making these approximations by solving 
the full Eqs. (f2"T)) - (f3"0")) for a stationary s. Results are 
given in the Appendix. One again finds that the thermal 
dephasing rate vanishes at zero temperature, and has the 
same high-temperature limit as Eq. (|4T[) . One also finds 
that corrections to the result Eq. (j4~4"|) are higher-order 
in 7. If we keep the ratio 7/A fixed and expand in 7, we 
find: 



K, th ] full = r Vtth + o( 7 3 ) 



(49) 



We now consider the additional effect of a driving force 
F(t) on the oscillator. The thermal factor Ai(t) consid- 
ered in the previous section will be unchanged by the 
presence of the force, assuming the initial state of the 
oscillator has variances consistent with being in equilib- 
rium. However, because of the drive, the means of the 
Gaussian (i.e. x and p) will be time-dependent; as a re- 
sult, the factor h.2(t) (c.f. Eqn. will also contribute 
to the suppression of Pti (*) an< ^ to the number splitting 
effect. To characterize this, we must calculate x(t) and 
p(t). 

We begin by writing: 



01 (t) 
a 2 (t) 



ip(t)) 
ip{t)) 



(51) 
(52) 



Note that in general, ai(t) ^ [a 2 (t)]* as x and p have 
complex parts; we only have ai(t) = [a 2 (£)]* at t = 0, 
when both x and p are real. Also note that with our 
rescaling, a± and 02 correspond to the usual a for a har- 



monic oscillator. As x 2 + p 2 — 4a2ai, we have: 



A 2 (t) = exp ( -2iX / dt'a 2 {t')ai{t') 



(53) 



Note the similarity to Eq. (5.16) in Ref. 2. In the positive 
P-function approach used in that work, it was important 
to have a different coherent states associated with the 
two states of the qubit (| |), | |)). In our approach, this 
"doubling" of degrees of freedom manifests itself in the 
fact that x and p have both real and imaginary parts 
(i.e. a 2 ^ («i)*). 

We use Eqs. (fl~5|) and (fl~6|) to obtain equations for a\ 
and a 2 - Assuming the drive is near the oscillator fre- 
quency fl, and that A <C 1, it is safe to make a ro- 
tating wave approximation to these equations. Writing 
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F(t) = Re [Fe luJdt ] , we thus find (recall we have set 

n = iy. 



dt al = 



dt a " 



—i iX 1 



l iX 

2 



2 

rr- + a 



oi + -F+(t)(54) 
02-^[F + (t)]\55) 



where F + (t) = Fe~' l ( UJd ~ 1 * >t . We see that due to the cou- 
pling A, the average variance cr x (t) + a p (t) acts as a time- 
dependent modulation of the oscillator frequency. This 
average variance is again completely determined from 
Eq. (|28|) of the previous section, and is completely in- 
dependent of the force F(t) or the initial value of a(t). 
In addition, as we have seen, it will have an imaginary 
part. 

The above equations are easily solved to yield: 



01 (t) 

02 (*) 



0i(t,O)ai(O) + - J dt' gi (t,t')F + (t') (56) 
OMO) -\J*dt'g 2 {tJ) [F + (t')}*(57) 



with 



<?i(M') 

92(t,t') 



e -i(t-t') e -y(t-t')/2 x 



exp 



-iA / dt" 
If 



a x (t") + a p (t") 



e *(t-t') e -7(t-t')/2 x 



(58) 



(59) 



exp 



- / A / dt" 



<7*(*")+*p(*") 



Eqs. (|53j) . (|56|) - (|59|) are enough to completely deter- 
mine the drive contribution Aa(t) to the decay of pn(t). 
In the limit of zero temperature, there is no thermal 
contribution to the decay of (i.e. there is no zero- 
temperature dephasing), and a x = a p = 1 for all times. 
This corresponds to the limit considered in Ref. [U; as 
we will see, the current approach yields the same answer 
found in that work. In the more general case of a non- 
zero temperature, the thermal number splitting will also 
influence the "drive" contribution to the number split- 
ting, through the dynamics of a x (t) + o~ p (i). 

To proceed, we make use of the small-7 solution for 
s(t) obtained in the previous section (c.f. Eq. (|52"j0. We 
stress again that this result would be exact if one started 
by making a rotating wave approximation in Eq. ([5]), 
as is commonly done in the quantum optics commu- 
nity. We also stress that corrections to this approxima- 
tion are negligible in the weak-damping, weak coupling 
(A«l,7<fl) limit of interest. We thus find: 



_ -i(l+Aa)(ti-fe) 



1 



Me -2i\at 3 

20^ (6°) 



»(l-Aa)(ti-t 2 ) 



1 - Me 
1 - Me~ 2iXat 



I -Me 



-2i\ati 



(61) 



The parameters M and a are temperature dependent, 
and are given by Eqs. (j3"3")) and respectively. Note 
that 31 and gi are not time-translational invariant, as 
they are sensitive to the initial time t = at which the 
coupling was turned on. Also note that in the zero tem- 
perature limit, we have M — ► and Xa — ► A — i"f/2. 
The propagators g\ and g2 then reduce to the expected 
propagators for a damped oscillator having a frequency 
of £1 ± A. The propagators also have this form at long 
times (ti7, *27 ^> 1) in the finite-temperature case. 
We can now write a\ (t) and a 2 (t) as: 



ai(t) = a 1 



02(f) 



a-2,s ■ e 



iu d t 



Oi(t) 
■02(f) 



(62) 
(63) 



Here, d\(t) and a 2 (t) describe the transient behaviour of 
ai and a 2 . This decomposition then allows us to write 
the drive-dependent contribution to Pn(t), A 2 (i), as: 

A 2 (t) = exp [-iA dr t - T Vydr t] x exp [Q tmns (t)] 

where 



i A dr — r Vi dr = — 2i\ai. s a,2,s 

Qtrans(t) = -2iX / dt' [ai (t')oa (*') 



(64) 
(65) 



+a M a 2 (i')e^ d '' + a 2 , s a 1 (i')e + ' i " d *' 

We thus see that the long time behaviour of ai, a 2 will set 
the drive contribution T Vt dr to the dephasing rate of the 
qubit (i.e. the exponential decay of pn at long times), 
while the transient behaviour of these two functions will 
be (presumably) responsible for the number splitting ef- 
fect. 



A. Results for the drive contribution to the 
dephasing rate 

In the presence of a drive, the total dephasing rate of 
the qubit (i.e. the exponential decay rate of p\i(t) at long 
times) will be given by: 



<p,th 



<p,dr 



(66) 



The purely thermal contribution T Vtt h is given in 
Eq. (|4"4"]l . Using the results of this section, we can cal- 
culate the drive contribution to the dephasing rate. One 
first finds the stationary amplitudes ai jS and a,2, s - Let- 
ting S = — ujd, we have: 



ai, s = 



02,. 



iF/2 
iS + iXa 

-iF*/2 
—i5 + iXa 



iF/2 



i(6 



A th ) + I 



■ tp,th 



-iF*/2 



-i(5- 



A th ) + I 



r 



ip,th 



(67) 
(68) 



Here, A t ^ is the thermal frequency shift of the qubit (c.f. 
Eq. (|4"3"]) h to lowest order in 7, A t ^ = A. Note that the 
stationary amplitudes are temperature-dependent! 
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Using these expressions, and letting j t h = (7+2r v , it / l ), 
we find that the drive-dependent contribution to the 
qubit dephasing rate takes the form: 



XAthith 



(S + A th f + {lfY\ [(5 - A th f + (Y) 

(69) 



Note that at zero temperature, A t h = A and T Vtt h = 0. 
In this limit, Eq. (|69|) reduces to the expression found 
m Ref. 0. We see that the effect of finite temperatures 
is to both provide an extra term in the dephasing rate 
of the qubit (i.e. T^th in Eq. (J44J)) , and also to mod- 
ify the form of the drive-dependent contribution to the 
dephasing rate. 



B. Results for the drive contribution to the 
dephasing spectrum 

As discussed, the transient parts of a\{t) and a 2 (t) will 
lead to oscillating contributions to p+\. Using the above 
results, one finds that these transients are given by: 



Si (t) 

02 (*) 

with 
Ai = 

A 2 = 

B x = 
B 2 -- 



1 



1 - Me- 2lXat 
1 



I -Me 



— 2iXat 



(70) 

A 2 ■ e*^ 1- Ao ^* + B 2 e i ^ d ~ 2Xa " ,t 

(71) 



(1 - M) x 
o(0) 



iF/2 



iS + iXa 
(1 - M) x 

a(o)r- ~ iF72 

-2iAaAf 
-2iAaM 



1 - 



2M 



iXa 



(1 - M) iS - iAa 



(72) 



(73) 



-i<5 + iAa 
iF/2 



2M 



iXa 



(1 - M) iS + iXa 



(iS + iXa) (iS — iXa) 

-%F*/2 
(iS + iXa) [id — iXa) 



(74) 
(75) 



The contribution to the number splitting from the 
transients of a\ and a 2 is given by Eq. (|65p . In gen- 
eral, the result will be a complex pattern of peaks, with 
no simple relation to the original photon number distri- 
bution of the mode. There are however several important 
situations where things simplify, as we now show. 



1. Zero temperature limit 

In the zero temperature limit, we have Bi = B 2 = 0, 
and: 



Ai 
A 2 



o(0) - oi, s 

HQ)]* - a 2. 



(76) 
(77) 



Only three terms are present in Qtrans 

AiA 2 



exp [Q transit)] = exp 

2a\.A 2 X 



1 







— 2i\at 



exp 
exp 



Xa — 5 
2d2, s AiA 
Xa + 5 



g*5t g— iXat 



2X 



' g—i8t g — iXat 



(78) 



Each factor above would yield a Poisson distribution 
of oscillating exponentials. The combination of all three 
would thus yield a complex distribution of peaks in the 
spectrum of p^, with weights being given as a convo- 
lution of three Poisson distributions. Ref. d considered 
an initial condition where the oscillator is in a state re- 
laxed to the qubit being down: a(0) = (o2,s)*- This is 
the correct initial state for an absorption experiment (c.f. 
Eq. 10) j where the coupling is always on, and the qubit 
is initially in its ground state. In this case A 2 = 0, and 
only the third factor in Eq. (f78|) contributes. One thus 
finds that Fourier transform of p-\i(t) consists of a set 
of Poisson-distributed peaks whose frequency spacing is 
approximately (in the small-7 limit) Aw = S + X. Note 
that if we assumed the spectrum of pn directly reflected 
the number statistics of the oscillator ( c.f. Eq. ^ ), we 
would expect a spacing between peaks of Auj = 2X. Also 
note that though one does obtain a Poisson distribution 
of peaks, the mean of this distribution does not corre- 
spond to |a(0)| 2 (i.e. the mean of the original number 
distribution before the measurement was turned on). It 
is unclear to what extent we can say we have "measured" 
the initial number distribution of the oscillator in this 
case. Note that if we specialize further, and chose the 
detuning of the driving force to be S = A, then the sim- 
ple expectation of Eq. ([7]) does hold : one obtains peaks 
with the expected spacing 2A, and whose weights corre- 
spond to a Poisson distribution with mean |a(0)| 2 . 

Let us now consider a different approach. Returning 
to the full expression, we see that only the first factor in 
Eq. (|75|) yields a set of peaks with the expected spacing 
of Aw = 2A, regardless of the detuning S. Consider the 
"strong measurement regime" , where A is much larger 
than both S and 7/2. Further, take the initial state of the 
oscillator to correspond to A = (i.e. a driven damped 
oscillator, but no coupling to the qubit). One then has 
|o(0)| 3> |o2,s|- In this limit, one can approximate 

Ai ~ a(0), A 2 ~ o(0)*. Only the first factor in Eq. (JTSJ) 
will contribute to the peak spacing; the second and third 
terms will be suppressed by the small factors |ai )S /a(0)| 
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and |ci2,s/a(0)| respectively. Thus: 



exp [Q transit)] ^ exp 



l«(0)| 2 



1 



?7 
2A 



(79) 



The Fourier transform of puit) will thus consist of a 
Poisson-distributed set of peaks with the correct spac- 
ing of Auj = 2A; moreover, in the small damping limit 
7 <C A, the mean of this distribution agrees with the 
mean of the original state, |ao| 2 . We have thus shown 
that at zero temperature, one can reliably extract the 
initial number distribution of the driven oscillator (arbi- 
trary detuning 5) from the time-dependence of pn us- 
ing the "number splitting" interpretation of Eq. . To 
do so, one must keep the qubit-oscillator coupling A off 
before t = 0; further, for t > 0, the coupling must be 
suitably large: A > S, 7. Note that one may still have 
Ac 1. It is worth emphasizing that the initial condition 
we consider here (i.e. Hi nt = before t = 0) is not appro- 
priate for the absorption experiment of Ref. 0- However, 
it could conceivably be realized in experiments involv- 
ing a nano-mechanical oscillator capacitively coupled to 
a qubit, as in such systems the qubit-oscillator coupling 
is proportional to a tunable voltage. 



2. Drive number splitting, finite temperature 

We now consider the transient drive contribution in the 
presence of a finite temperature. We again consider the 
"strong measurement limit" A > (5, 7. In this limit, the 
initial amplitude of the oscillator |a(0)| is much larger 
than |ai jS | or |o2,«|< Thus, the only significant term in 
Q transit) is the one proportional to A1A2, as this is the 
only term which scales as |a(0)| 2 . Using A1A2 ~ (1 — 
M) 2 |a(0)| 2 (c.f. Eqs.JT^-GSD), we find: 



t e -2i\at' 



Qtrans = ~2i\AlA 2 dt' 

Jo (I - Me 2 

K0)l 2 



-2i\at 



1 



a 



y 

— 2i\at 



(80) 



1 + M- 



1 



1 



Me -2z\at 

(81) 



If we now combine the above results for the drive con- 
tribution to Pi lit) with the thermal contribution Ai(t) 
(c.f. Eq. (36])) we find: 



p u it) = e -i^t e -i^4rt e -r^ dr t x 

Pik = -2\a;T eff ,F eff ,5) 



(82) 



form of this distribution is given in the Appendix. To see 
the correspondence most directly, make the substitutions 
e ~/3f2 _^ j^r anc j _2Aai in the generating function 
given in Eqs. (|B15|) - ([B17]) . We find that in the small 
damping limit, T e // and P e // appearing in Eq. (|83[) are 
almost identical to the initial temperature and driving 
force on the oscillator: T e ff is again given by Eq. ([40]) 
and 



F 



F 



ff 



F 



8n 



cqj 



(84) 



Note the close connection between Eq. ([83]) and the 
naive "number splitting" expectation for p\iit) given in 
Eq. ([7]). We have thus shown that in the "strong mea- 
surement limit" A > 7,i, the spectrum of p\i(t) accu- 
rately reflects the number statistics of the mode: in this 
limit, each peak in the spectrum of p^ corresponds to a 
number state in the mode, and the corresponding peak 
weight to Pin). The width 7„ of the peak corresponding 
the photon number n is easily found to be: 



In = T Vt dr + (2n + 1) (r v ,th + ^) 



7 
2 



(85) 



In Fig. 1, we show results for the dephasing spectrum 
Siu>) for various values of the qubit - oscillator coupling 
strength A, where: 



S(u) 



dte i»t(pnif) e ixt 
W(o) 



(86) 



Note we have shifted the spectrum so that the peaks as- 
sociated with number states should appear at frequencies 
u> equal to even multiples of 2A. One also expects from 
Eq. that the weights of the different peaks peaks in 
the spectrum should correspond to the oscillator num- 
ber distribution. These weights may be approximately 
obtained from 5(tj) via: 



P(m) 



f(2m+l)A 



approx 



1m J(2m-1)A 



dujSiu) (87) 



where h m is the maximum of Siu>) in the interval 
[(2m — 1)A, (2m + 1)A)]. Number distributions obtained 
from Siu>) using this procedure are shown in Fig. 2. 
Note that one could be more sophisticated in extracting 
the weights of peaks from we use Eq. (|57|) only for 

simplicity. 



Y / Pin;T eff ,F effl S)e- 2 ^i n+ i) t (83) 



Here, P(n;T, F,5) is the photon number distribution of 
an oscillator at temperature T which is driven by a force 
of magnitude F and detuning S; P(fc) = Y^ n e%k ^'i n ) 
is the corresponding generating function. The explicit 



VI. CONCLUSIONS 

In this paper, we have used a conceptually simple 
Wigner function approach to study the physics of a qubit 
dispersively coupled to a driven, thermal harmonic mode. 
In general, back-action effects prevent the qubit 's dephas- 
ing spectrum (i.e. the Fourier spectrum of Puit)) from 
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FIG. 1: Fourier transform of the dephasing spectrum, S(u>) 
(c.f. Eq. (JHU) for different values of A. For each plot, we 
have taken the oscillator damping 7 = 10 -4 , a driving force 
of magnitude \F\ = 37 and detuning 8 — 0, a temperature 
which corresponds to n eq = 1, and an initial condition which 
corresponds to the driven thermal oscillator at zero qubit- 
oscillator coupling. Based on the number splitting interpre- 
tation of Eq. 0, one expects to see peaks in S(u>) whenever 
uj = 2A. For small values of A, S(u>) is in poor correspon- 
dence to the number distribution; note that there are addi- 
tional peaks occurring roughly at u) equal to odd multiples of 
A. These unwanted peaks become increasingly insignificant 
as A is increased. The weights of the peaks in plot c) are in 
good agreement with the expected number state distribution 
for the initial oscillator state (see Fig. 2). 
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FIG. 2: Number distributions derived from the dephasing 
spectrum, using the same oscillator parameters as in Fig. 1. 
The black bars correspond to the actual photon number dis- 
tribution associated with the initial oscillator state. The re- 
maining bars correspond to distributions obtained from the 
dephasing spectrum via Eq. (|87|l . For A = 0.01, there is a 
good correspondence. Note that the deviations from the true 
distribution in this case are mainly due to the inaccuracy of 
using Eq. (|87[1 to obtain the weight of peaks. 



APPENDIX A: CORRECTIONS TO THE 
THERMAL DEPHASING RATE 



In the main text, the thermal dephasing rate is cal- 
culated starting from the master equation Eq. ([9]), using 
a weak-damping approximation. As discussed, this ap- 
proximation is equivalent to having dropped terms pro- 
portional to a 2 and (c^) 2 in Eq. ©, as is often done in 
the quantum optics literature; this approximation was 
also used in Refs. Q], 0- In this section, we discuss correc- 
tions arising when one goes beyond this approximation. 



1. Results from the full, non-Linblad master 
equation 



being a perfect mirror of the number statistics of the ini- 
tial oscillator state. We have shown however that for 
a sufficiently strong qubit-oscillator coupling, one can 
indeed reliably extract the photon distribution of the 
driven, thermal mode from the dephasing spectrum. We 
have also calculated the oscillator-induced dephasing of 
the qubit in the case of a thermal, driven mode. 

We thank Alexandre Blais and Mark Dykman for use- 
ful conversations. This work was supported by NSERC, 
FQRNT and the Canadian Institute for Advanced Re- 
search. 



If one retains the counter-rotating terms in Eq. ([5]), 
one finds corrections to the result Eq. ([44]) . These correc- 
tions are irrelevant in the small damping, small coupling 
(7, A< 1) limit of interest. Explicitly, expanding T^^h 
in 7 for a fixed A, one finds that corrections to Eq. ([44]) 
only arise at order j 3 : 



\p<p,th]f u u 



in eq - 

7 3 1-3A 2 +4A 4 
4A^ (1-A 2 ) 2 

+0( 7 5 ) 



(Al) 

n eq (l + n eq )(l + 2n eq ) 
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Comparing against the corresponding expansion of the 
approximate answer, Eq. (|46[) . we see that the correction 
terms arising from using the full equation are suppressed 
by a factor A 2 compared to the order-7 3 terms in approx- 
imate answer.. 

Similarly, if we expand the full answer to lowest order 

A, we find that corrections to Eq. l[4"8[) are suppressed by 
2. 
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\F<e,th]f ull - A 



An eq (l + n eq ) 



-)+0(X^ 



(A2) 



Results for a counter-rotating qubit-oscillator 
interaction 



One can also calculate the thermal dephasing rate of 
the qubit using an qubit-oscillator interaction Hamilto- 
nian which retains counter-rotating terms, c.f. Eq. ([50)) . 
It was recently sug gested that such terms could play a 
significant role [151 ]. Repeating the treatment which led 
to the full dephasing rate above (i.e. no small damping 
approximation), we obtain: 



[T v ,th] H - 



-Re 

2 




HX 



'eq 



-r_i 



where 
.4 



^2 - 7 2 + 2y/l + Ai 1 X{2n eq + 1) - 4A ; 



(A3) 



= (A4) 



Note that this expression is identical to the previous 
Eq. ([4"4"[) except for terms which are small as 7 or A. To 
see this explicitly, we can expand this result for small A 
to obtain: 

[T v ,th} H . nt2 = 
A 2 



7 



4n eq (l + n eq ) + 7 (1 + 2n eg ) z + 0(A 4 )(A5) 



The second term above corresponds to the correction 
coming from using the non-rotating wave interaction 
Hamiltonian; it is suppressed by a factor 7 2 . Thus, as 
claimed, the non-rotating terms in the qubit-oscillator 
interaction Hamiltonian Hi nt play no role in the weak 
damping, weak coupling (A <C 1,7 <C O) limit of in- 
terest. Note that T Vtt h for an interaction Hamiltonain 
H-int^i was also calculated in Ref. fl5l : they found a result 
similar to Eq. I|A5[) above. 



APPENDIX B: NUMBER STATISTICS FOR A 
DRIVEN THERMAL MODE 

In this appendix, we calculate the generating function 
for the number distribution for a driven harmonic oscil- 
lator which is also coupled to an equilibrium dissipative 
bath. The Wigner function for such a state will be a 
Gaussian, with the variances determined by equiparti- 
tion. The means of the Gaussian will follow the equations 
of motion for a damped, driven oscillator. In general, we 
can represent this state in the form: 



Pdr 



= DpoD- 1 



(Bl) 



where po is the equilibrium density matrix of the oscilla- 
tor, and D is the displacement operator: 



D = exp (a^a) — a^a 



(B2) 



Finally, ao is the coherent state amplitude corresponding 
to the driven, damped motion of the oscillator; only its 
magnitude will be important for the number statistics. 
We have: 



"ol 



|F| 2 /4 



(^-fi) 2 +7 2 /4 



(B3) 



To proceed, it is useful to represent po using a P func- 
tion: 



Pa = J daP (a)\a){a\ 
1 

7T 



(B4) 



P (a) = -( e ^ a -l)exp(-|a|V a -l)) (B5) 

(B6) 



It thus follows: 



p = J daPo(a)\a + ao)(a + ao\ (B7) 

The number distribution is then given by: 

P(n) = (n\p\n) (B8) 

|a+ an | H« + «0| 2 " (B9) 



= J daPo(a)e~ 



00 plix 

rdr I d6 
11 ^0 



r 2 e^ n -2r|a | cosfl- 



^2 (r 2 + |a | 2 + 2ra cos6>)" 
(BIO) 

The corresponding generating function is thus: 
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P(k) 



= j2 p ( n V kn ( BU ) 

n 

f J3u> _ 1 \ Z" 00 /•2tt 

= exp[-|a | 2 (l-e lfc )] ^ J y rdr J d8 cxp [-r 2 (e pu - e ik )] exp [-2r|a | cob 0(1 - e lk )] 



f-OO 

cxp[-|a | 2 (l-e lfc )] 2(e /3w -l) / rdr cxp [-r 2 (e^ - e 4fe )] J [2r\a \(l - e lk )\ 
I Jo 



The first exponential factor above corresponds a gener- 
ating function for a Poisson distribution; this would cor- 
respond to the number statistics of a driven, zero tem- 
perature oscillator (i.e. coherent state statistics). The 
second term corresponds to the effect of a non-zero tem- 
perature. In the absence of the Bessel-function factor, it 
would yield the generating function for a thermal distri- 
bution. In the presence of a drive, the Bessel-function 
factors modifies the form of this term. Thus, not surpris- 
ingly, the number statistics of a driven thermal oscillator 
is not simply a convolution of a Bose-Einstein distribu- 
tion with a Poisson distribution. 

We can evaluate the expression for P(k) further by 
making use of: 



> z 2l 



(B12) 



and 



H drr 1+2m e- Ar2 = —A^ 1 (B13) 
Jo 2 



Some algebra thus yields: 



P(k) = Ptherm{k) X P drive 

(k) x 

Pcorr 

(k) (BU) 



with 



Ptherm (&) 

P d (k) 

Pcorr(k) 



1 - er fm 

I _ e -pn e ik ( B15 ) 
exp[-|ao| 2 (l-e jfc )] (B16) 

|Qo1 6 ,(l-e* k ) 21 
t£ l\ (1 - e-We^) 1 

(1-e 



E 



= exp 



\a \ 2 e~^ 



(B17) 
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